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PREFACE 


A Float Zone Workshop was held at MSFC on September 22-23, 1981. The pur- 
pose of the workshop was to bring together key industrial, academic, and government 
scientists to (1) Determine the possible advantages from space processing (by float- 
zoning) and where these improvements are needed by industry and government and 
(2) Forecast how much power is required for float- zone space processing. 

The workshop was structured into two days. The first day consisted of reports 
from various industrial and government producers and users of float- zone material. 

The second day was the presentation of present research efforts in the float- zone 
area. 


The list of attendees and an agenda are given in Tables 1 and 2. 
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INTRODUCTION 


E . L . Kern 
Chairman, FZWG 
Del Mar , CA 


The purpose of the September 1981 Float Zone Workshop was to bring together 
key industrial, academic, and government scientists to address the questions of: 

1) What are the needed improvements in solid crystals? 

2) What is the advantage of processing in space? 

3) How much power will be needed? 

Based upon this rationale, the present (and soon to be) principal investigators 
presented their recent results and future plans for providing the necessary ground 
based research. 

The set of presentations which addressed the above questions was preceded by 
a current status of the Material Processing in Space (MPS) program by Dr. Bill Oran, 
of that program office at NASA Headquarters. A review of past space efforts was 
provided by a chronology of the types of space vehicles and experiments that were 
carried out, presented by Gene McKannan of the MPS Proiect office. Dr. Robert 
Naumann and Professor Bill Wilcox then reviewed what we have learned from past 
experiments . 

My introductory figures, as Chairman of the Float Zone Working Group (FZWG), 
follow. Figure I-l presents my definition and purpose of the FZWG; Figure 1-2 
defines the purpose of this particular workshop. The need to forecast power needs 
in the future was outlined in Figure 1-3. Past, present, and proposed contractual 
or grant programs are listed in Figures I- 4a and I- 4b, with several of these providing 
results in the papers included in this report. I recently returned from reviewing 
floar-zone activities in Europe, by attending the 32nd Congress of the International 
Astronautical Federation and visiting two of the investigators in Germany (Drs. 

Kolker and Eyre). The efforts are listed in Figures I-6a and I-6b. Figure 1-7 gives 
a proposed program roadmap for the next several years, as a means of defining when 
our ground-based investigations should be able to lead to a definition of major flight 
experimentation and hardware. 

My impression of the European efforts is that there are several good research 
programs defining the limits of bouyancy convection and Marangoni flows. These 
studies are not being integrated into the silicon float -zone experiment of Spacelab 1 
and thus a solid ground-based foundation is not being provided for this experiment. 

In particular, a high thermal gradient will tend to maximize Marangoni flow and the 
me It -solid system thereby cannot be characterized to the extent of indicating the 
advantages of space zoning of silicon. Although we will not be zoning silicon in 
space as soon as Spacelab 1, our ground modeling and characterization will lead to 
more definitive experiments . We will remain in communication to learn of any of 
their results which would apply to our studies. 
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Figure I-l. 
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Figure 1-2. 
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Figure 








A. NASA: MATERIALS PROCESSING IN SPACE OVERVIEW 


William Oran 
NASA Headquarters 
Washington, DC 


Dr. William Oran, from the Materials Processing in Space (MPS) Program Office 
at NASA Headquarters, welcomed the workshop attendees. He presented a general 
overview of the MPS program , as outlined in Table I-l. Because of the present 
budget reviews and possible reductions in the MPS budget. Dr. Oran indicated that 
the size and speed of progress on the float- zone program could not be predicted at 
this time (September 1981). 
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TABLE I-l. MATERIALS PROCESSING IN SPACE 


Research Area 

Research Task 

^ 

Space Hardware 

Current Program in Physical and Chemical Sciences 

Containerless Science 

Glass Fining Studies 

Acoustic Containerless 
Experiments System 


Positioning Force Interaction 
Studies 

Acoustic Containerless 
System 



Electromagnetic Container- 
less Experiments System 

Fluid Behavior 

Crystal Growth Kinetic Studies 
and Fluid Dynamics Influences 

Fluid Experiments System 


Low Gravity Fluid Behavior 


Combustion Science 

Gravity Influences on Flame 
Initiation, Spreading, Kinetics, 
and Two-Phase Combustion 

Planned 

Cloud Physics 

Cloud Formation Processes, 
Ice /Water Role, Charging 
Phenomena, Chemistry 

Planned 

Vacuum Science 

Characterization of Gas 
Molecular Density and 
Distribution Behind Orbiting 
Shield Configurations 

Planned 

Current Program in Materials Processing 

Crystal Growth 

Semiconductor Crystals 
Grown from the Melt with 
Uniform Compositions 

Solidification Experiment 
System 


Semiconductor Crystals 
Grown from the Vapor with 
Controlled Morphologies 
and Low Defect Densities 

Vapor Crystal Growth 
Furnace 

Solidification 

Metal Alloys with Aligned 
Phase Structures 

Solidification Experiment 
System 


Composite Variable Density 
Solids with Uniform Distri- 
bution of Second Phase 

Solidification Experiment 
System 

Containerless Processing 

Glass Shell Formation and 
Uniformity Studies 

Acoustic Containerless 
Experiments System and 
Drop Dynamics Module 


Glass Formation in New 
Materials Systems 

Acoustic Containerless 
Experiments Systems 

Chemical Processing 

Formation of Monodisperse 
Polymer Latex Spheres Ki- 
netics , and Two- Phase 
Combustion 

Polymer Latex Reactor 
and Experiment Kinetic 
Chambers 


12 




B. REVIEW OF PAST FLIGHT EXPERIMENTS PERTINENT TO FLOAT-ZONE 


Gene McKannan 
NASA-MSFC 
Hixntsville, AL 


Mr. Gene McKannan from the Materials Processing in Space Project Office pro- 
vided a review of relevant material processing experiments that were done in space 
since 1973. These included experiments carried out in Skylab, in the Apollo- Soyuz 
mission, and on SPAR sounding rockets. The lists of experiments discussed are 
shown in Tables 1-2 and 1-3. 
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TABLE 1-2. MATERIALS PROCESSING IN SPACE 
1973 FLIGHT EXPERIMENTS PERTINENT TO FLOAT-ZONE 
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TABLE 1-3. MATERIALS PROCESSING IN SPACE 
1975 FLIGHT EXPERIMENTS PERTINENT TO FLOAT -ZONE 
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C. THE NATURE OF MELT PROCESS CHANGES IN MICROGRAVITY 


William Wilcox 

Dept . of Chemical Engineering 
Clarkson College 
Potsdam, NY 


The lack of bouyancy- driven convection in melts and the gas phase micro grav- 
ity will have pronounced effects on impurities within the melts, such as the case for 
doped silicon crystals. The change in growth characteristics could be complicated 
and significant and are difficult to predict. The first change is that a very thin 
adsorbed impurity layer could retard convective flow because of surface tension 
gradients (Marangoni flow) , whereas the addition of bulk melt convection at g = 1 would 
not allow such easy flow retardation. 

The second effect is the lack of convective flow on the transport of volatile 
impurities to the melt surface, thus reducing evaporation as compared to g = 1 proc- 
essing. Normally volatile impurities (such as As or In in silicon) need much lower 
overpressures to prevent evaporative losses. 

The third effect is the lack of convective flow in the gas next to the melt. 

This will change the dependence of evaporation rates on pressure and could change 
the heat loss magnitudes significantly. 

These process changes can only be observed in space experiments and point to 
the need for precusory space experiments to be made soon. 



Thermocapillary Convection 

Surface film or layer - stops all convection. 

Impurity adsorption (monolayer or less) 

Retards flow increasingly as: 

1. Solubilities in bulk phases decreases. 

2. Amount of impurity increases. 

3. Impurity dependence of surface tension increases. 

How to observe: 

1. Rate of bubble rise decreased. 

2. Do not have rapid movement of surface, as determined by flowing 
particles on surface. 

REDUCED EVAPORATION 

No convective transport of volatile species. Diffusion only. Consider 1 
dimensional transport through inert gas layer of thickness X at total pressure P 
condensing on surface so cold that vapor pressure is zero. Vapor pressure p at 
melt surface. Evaporation flxix is: 

N = - LN (1 - £) MOUM^S 


If zone contains n moles and has free surface area A, then in time t fraction 
evaporating is: 


F 


DPAt 

RTXn 


LN (1 
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EXAMPLE 


2 

DP = 0.1 cm atm/s 

Zone: 2 cm diameter, 2 cm long, containing 0.5 mol 

T = 1600 k 
X = 1 cm 

p/p Evaporation 

0.5 5%7hr 

0. 1 0. 7%/hr 

0.01 0.07%7hr 

0.001 0.007%7hr 

To hold losses below 0.1%, need overpressure about 100 x larger than vapor pressure. 
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On Earth, evaporation rates are much larger and much less pressure dependent — 
more like (laminer) or (turbulent). 

May have condensation in vapor to form mist or fog. 

Would increase evaporation rate and reduce radiation losses from zone. Re-evaporation. 
Complicated 
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II. CRUCIAL MATERIAL NEEDS FOR 
DEVICES AND THE ROLE OF SPACE 




A. THE NEED FOR GROWING CRYSTALS IN SPACE 

E. L. Kern 
Constiltant 
Del Mar, CA 


I have summarized some of the areas of growth understanding and some of the 
payoffs to microgravity crystal growth in Figures II- 1, II- 2, and II- 3. Comparison 
of the real zone shape at g = 1 (Fig. II-4) to the ideal shapes possible in microgravity 
(Fig. II- 5) shows the flexibility to be gained by processing in space. The possible 
heating methods are reviewed (Fig. II- 6). Figure II- 7 lists critical devices which 
need the uniformity and lower defect density which might be obtained by growth in 
space. Another payoff is to grow crystals which are difficult or impossible to grow 
on Earth, such as Silicon- germanium alloys and III-V and II -VI compounds (Figs. 

II -8 and II- 9). 
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Figiire II- 
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Figure II- 
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Figure II- 
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Figiire II- 
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Figure II- 
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B. MATERIAL AND PROCESSING NEEDS FOR SILICON 
SOLAR CELLS IN SPACE 


Henry W. Brandhorst, Jr, 
NASA -Lewis Research Center 
Cleveland , OH 


Dr. Brandhorst presented the NASA technical concerns on space grade solar 
cells. His first figure projects the rather low and constant solar power needs, until 
space platforms (PEP) increase the need to 50 kW new production /year in the later 
1980's. 

The main technical concerns are cell efficiency and radiation performance reduc- 
tion (Figs. II-ll and 11-12). The main causes for radiation degradation are complexes 
of impurities and crystal defects (Fig. 11-13). Vacancy -oxygen- carbon complexes 
form the largest concentration of centers which dissapated carriers within the solar 
cells. Boron -oxygen vacancy complexes and Boron-oxygen complexes also degrade 
the cell performance. 

Float- zone silicon is superior to the normally used silicon made by the 
Czochralski process. Figure 11-14 shows a direct correlation, with the very low 
oxygen leading to the improved float- zone silicon performance. Figure 11-15 shows 
that reducing oxygen and carbon (as float- zoning does) allows for much faster low 
temperature annealing of the defects. Figure 11-16 shows a marked improvement 
using gallium as the dopant in float -zone, as opposed to boron. 

The advantages of improved crystal purity and perfection upon cell properties 
are summarized in Figure 11-17. 
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CUMULATIVE PAST AND PROJECTED 
SOLAR CELL POWER FOR NASA PROGRAMS 



Figure II -10. 
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RADIATION INDUCED DEFECTS IN P-TYPE SILICON 
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Figure II-IS. 



EFFECT OF MATERIAL ON SILICON SOLAR CELL PERFORMANCE 



1-MeV ELECTRON FLUENCE (E/cm^) 






TEMPERATURE (° C) 




EFFECT OF DOPANTS ON CELL DEGRADATION AND ANNEALING 



ISOCHRONAL ANNEALING TEMP. (° C) 



SILICON MATERIAL AND PROCESSING REQUIREMENTS FOR INCREASED 



IRRADIATION (POTENTIAL TO ATTAIN < 200° C). 



C. FLOAT-ZONE SPACE PROCESSING RATIONALE 


Lubek Jastrzebski 
RCA 

Princeton, NJ 


Because float- zone growth in space is considered for commercial applications, 
there is a tendency to first establish which devices or circuits will feasibly be made 
from space grown float- zone silicon from an economical point of view. I believe that 
this is not totally justified. The majority of material problems which affect the 
performance of the circuits can presently be eliminated by adjustment of processing 
or changes in circuit design. The best example is striations observed in CCD 
imagers , which one and a half years ago were considered to be an insurmountable 
problem caused by materials. They have recently been eliminated by changes of 
processing and circuit design. This indicates that any material problems which right 
now seems to be an obstacle could be solved by future developments in processing 
or circuit design in the period of a few years, and it will take this long before float- 
zone crystals of at least 3-in. diameter will be grown in space. Therefore, focusing 
of research effort in space right now on any particular problem could be a little 
premature. 

It seems that the first objective of the present program should be to demon- 
strate that growth in space will result in an improvement of silicon crystal properties. 
It is difficult to predict requirements for silicon of the future, because they depend 
on the directions in which technology will evolve and only some guesses on this 
subject can be made. It seems that future ICs will not need to interface with TTL 
logic, and therefore high resistivity material will be used for processing. In this 
material, the problem of oxygen thermal donors, and especially their inhomogeneous 
distribution, could become a critical issue (Fig. 11-19). It has to be pointed out 
that this problem cannot be simply solved by elimination of oxygen from silicon 
crystals because oxygen in silicon also plays a beneficial role (Fig. 11-18). In this 
case float- zone silicon crystals doped with oxygen and grown in space can be of 
interest . 
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THE ROLE OF FLOAT -ZONED SILICON IN ICs 


Material (Silicon) for IC's 


CZ - Majority 


( + ) 

Oxygen pinning (mechanical strength) 
- Oxygen gettering (decontamination) 


( - ) 

Defects (crystal growth) 
Low resistivity - problem 


Large change of resistivity 
possible. 


FZ - Special Applications - Easy dislocations 
(CCD imagers) movement 

- Lack of gettering 

- Inhomogeneity 


- No crystal growth related defects 

- High resistivity 

- Small changes of resistivity 


16 3 

Present silicon for IC's: Resistivity, few 0cm ('^ 10 cm" ) 

reason: circuits are designed for 5V to interface with TTL logic, 

future: computer on a chip (e.g. , 32 bit microproc.) interfacing with TTL 

logic will no longer be required — > lower voltages (speed) > 

higher resistivity material 100-1000 0cm range. 


Possible problem: Nonuniform oxygen distribution. 


Figure 11-18. 



CRITICALITY OF OXYGEN AND DOPANT RADIAL 
PROFILES IN MOS DEVICES 


Present Future (?) 

— — — level required 

for devices 

impurities 

oxygen donors 

'W\A 





CMOS 


CMOS, latch up, threshold voltage 
NMOS, dynamic, threshold control 



\ Pi 


tou OK. do 

CouuXc 



Present: inhomogeneities can be overcome - inconvenience. 


ox, 

oU)UOVS 

CoxAttulv<iiwti 


Future: possible problem with oxygen, uniformity — problem of point defects 

and their nonequilibrium trapping (their effect on initial stage of 
oxygen precipitation kinetics (?)) 


Silicon requirements in the future — depends on direction in which system 
development will go. 

Flat Radial Profile of dopant required. 


For mass production space advantages: 

crystal growth - size of the grown crystals 
processing - wafer thickness 
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Figwe II-IR. 



FLOAT -ZONED SILICON FOR BULK DEVICES 


Float Zone Silicon - Power Devices - High Resistivity 
For very deep junctions - bulk silicon required. 

Nonuniformity of resistivity on micro scale - junction depth variations - 
breakdown - (due to electric field nonunifoimity) . 

Question: cost 

Power devices require a large amount of silicon material (compared with IC's) 
cost low (IC's) - material cost is a large contributor to device cost. 

Will improvement in device performance pay for additional cost? 

IR Detectors ^ ^ supporting IC's speed. 


Figure 11-20. 



D. GROWTH OF III-V AND II-VI ALLOY CRYSTALS 
IN SPACE BY VERTICAL ZONE LEVELING 

Klaus J. Bachmann 
North Carolina State University 
Raleigh, North Carolina 


Because of the change in density of zincblende structure compounds, e.g. , 

InP , upon melting, a vertical zone melting process proceeds in outer space without 
contact of the molten zone with the container walls providing thus the advantages 
of float- zoning and confinement of high decomposition pressures which cannot be 
obtained under the conditions of ground based research (compare Figs. 11-21 through 
11-28). This advantage is expected to result in substantial improvements in the 
perfection of important electronic materials, e.g., CdTe and InP, when grown by 
vertical zone leveling in space and promises, in addition, to provide unique and 
essential conditions for the growth of III-V and II-VI alloy crystals, e.g., 


Cd Hg^ Te and InP As^ 

X °l-x y 1-y 


An evaluation of various ternary and quaternary III-V 


alloy systems in terms of their suitability for experiments of very limited duration. 


and we suggest some work on InP As, and Cd Hg, Te that are of interest in the 

y 1-y X ^1-x 


context of near infared devices for optical communications and of far infrared 


detectors for mosaic focal plane sensors, respectively. For details we refer to a 


recent review: K. J. Bachmann, H. Schreiber, Jr., and F. A. Thiel, Progr. in 

Crystal Growth and Characterization, Vol. 2, pp. 171-206 (1979). 
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Figure II - 28 . 



E. AIR FORCE FUTURE TECHNICAL REQUIREMENTS 

FOR SILICON 


Patrick M. Hemenger 
Wright -Patterson Air Force Base 
Wright -Patterson AFB , OH 


The Materials Laboratory (ML) is presently supporting materials improvement 
programs that include a range of semiconductors (e.g. , Si, GaAs, HgCdTe, InAs) 
in both bulk and thin layer forms which are aimed at solving a spectrum of problems 
ranging from impurity control to radiation hardening of integrated circuits. This 
presentation emphasizes the Si based bulk materials because they are more directly 
related to this meeting. However, closely related work on other materials and on 
thin-layer Si are summarized for completeness. Figures 11-29, 11-30, and 11-31 
summarize the three topics covered ; 

1) Si based materials requirements 

2) ML supported programs 

3) ML inhouse characterization facility. 

Several general thoughts and questions come to mind when considering the 
growth of crystals in space : 

1) Some of the desired improvements for bulk Si should be realized directly, 
namely higher purity, greater uniformity, and better overall control of defects and 
dopants . 

2) Can the Si-Ge system (Si-rich) be grown with the quality necessary for 
use in wavelength tuned detectors, both intrinsic and extrinsic? 

3) Will magnetic field growth retain any advantages? 

4) Can the solution growth process be improved in the zero gravity environ- 
ment if it allows a larger solution zone while maintaining stability? 

The Air Force and the ML plan to participate actively in this project. This 
participation will include, on a continuing basis, characterization of the space grown 
crystals and evaluation of their applicability to Air Force systems. The complete 
inhouse characterization facility which is outlined on the attached sheet will be made 
available for this project. 
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REQ UIREMEN T S /B AC KGRO UND 


HIGHER QUALITY DOPED SILICON IS REQUIRED SUITABLE FOR 
FABRICATING MONOLITHIC FOCAL PLANE ARRAYS WITH ON CHIP CCD 
SIGNAL PROCESSING. 

-Gallium Doped Silicon with Greater Uniformity of 
Responsivity and Higher Responsivity is Required 
for the Next Generation Tactical FLIR. 

-Indium Doped Silicon with Higher Operating Temperature 
and Greater Uniformity of Responsivity and Higher 
Responsivity is Required for Very Large Focal Planes of 
the HALO Type for Space Surveillance. 

-Ultra-high Purity Silicon Required for Laser Guided 
Weapons Systems. 

-Very Long Wavelength Doped Silicon Material Such as 
Si: As is Required for Space Missions. 

OTHER SILICON BASED AIR FORCE MATERIALS REQUIREMENTS 

-Bulk Silicon with Larger Minority Carrier Lifetime and 
Thin Silicon Films ( 2000^) on Sapphire with Higher 
Mobilities Required for VLSI Applications. Air Force 
Program all Silicon Based. 

-Long Lifetime Si-based Microwave Sources and Detectors 
are Required for Phased-array Radars (Over the Horizon 
and Spaceborne) In Early Warning Systems. 


Figure 11-29 



Si-Ge ALLOY INFRARED DETECTORS 

OBJECTIVE 

• DEVELOP A DETECTOR MATERIAL SYSTEM TO EXACTLY "TUNE" 
THE SPECTRAL RESPONSE TO NECESSARY OPERATIONAL CON- 
DIT IONS 

APPROACH 

• MANUFACTURE IR DETECTORS FROM SILICON-RICH Si-Ge 
ALLOYS 

PAYOFF 

• FOR EXTRINSIC DETECTORS (FLIRS) A HIGHER OPERATING 
TEMPERATURE FOR ANY SPECIFIED ATMOSPHERIC WINDOW, 
THEREBY REDUCING POWER REQUIREMENTS, COOLER COST 
AND WEIGHT 

• FOR INTRINSIC DETECTORS (RANGEFINDERS), A FACTOR 
OF 2 INCREASED SENSITIVITY AT 1.06 MICRONS, WHICH 
PRODUCES AN EVEN LARGER FACTOR IN TARGET ACQUISITION 
RANGE 

Figure 11-30. 
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FACILITIES /CAPABILITIES 


• electrical and optical 

-HIGH RESISTIVITY HALL FACILITY 

-PHOTO-HALL FACILITY 

-LASER ANNEALING FACILITY 

-IR ABSORPTION SPECTROMETRY 

-PHOTOLUMINESCENCE FACILITY 

-PHOTOCONDUCT IV IT Y/PHOTOTHERMAL 
IONIZATION 

-THEORETICAL CAPABILITY 


• SURFACE /CHEMISTRY 

-SCANNING AUGER MICRO 
SCOPY 

-X-RAY PHOTOELECTRON 
SPECTROMETER 
-ION SCATTERING SPECT 
ROMETER 

-SECONDARY ION MASS 
SPECTROSCOPY 
-SCANNING ELECTRON 
MICROSCOPE 

-AUGER DEPTH PROFILE 
INSTRUMEI^T 
-ELLIPSOMETER 
-LASER-RAMON SPECTRO- 
METER 

-LOW ENERGY ELECTRON 
DIFFRACTION 


PLANNED 

-CDLTS & C-V MEASUREMENTS FOR FILMS 

-AC TECHNIQUE FOR CHARACTERIZING WAFERS/FILMS 

-UPGRADE OF EXISTING FACILITIES/CAPABILITIES 


Figure 11-31. 
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F. INDUSTRY NEEDS FOR SILICON CRYSTALS AND STANDARDS 


Kenneth E. Benson 
Bell Laboratories 
Allentown, PA 

(Chairman ASTM F-1 Committee on Electronics) 


The trend of the device fabrication industry requiring larger crystals is 
reviewed. The standard Czochralski wafer size is now 100 mm (4 in.) and going to 
125 mm (5 in.) within 2 or 3 years, with 5-in. wafer device lines now coming "on- 
stream" (Fig. 11-32). The mass of Czochralski crystals has steadily increased over 
the past years and crystals of 30 kgm mass are now standard (Fig. 11-33). 

A real need of the device industry, which cannot be met by the nonuniformity 
in production crystals, is for measurement standards. This is an area where the 
NASA crystal growth program could make an important contribution to industry and 
government needs. Standards are needed for resistivity (four-point probe and 
spreading resistance) and for chemical composition of oxygen and carbon impurities. 
The ranges of properties and uniformities for these standards are given in the 
following figures. 
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STANDARD REFERENCE SILICON MATERIALS NEEDS 


ASTM Committee F-1 on Electronics 

• Resistivity Standards 

Four -Point Probe 
Spreading Resistance 

• Oxygen Standards 

• Carbon Standards 


RESISTIVITY STANDARDS 

Nominal Doping Ranges — for N and P Type Material, 
<100> and <111> Orientations 5 x 10^^, 1 x 10^^, 

1 X 10^ , and 1 X 10^ Jcm^ 

Uniformity <0.5% 


OXYGEN STANDARDS 
Nominal Concentration — 0.1, 1, 40 ppm 
Uniformity ±5% 


CARBON STANDARDS 
Nominal Concentration — 0.1, 1, 10 ppm 


Uniformity ±5% 



125mm 



ooooo oooo 
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A) RANGE - 0.0008 TO 5,000 fi-cm 

B) VOLUME SAMPLED - 0.02 cm3 

C) PRECISION ±1% 




SPREADING RESISTANCE PROBE 

WAFER AND IN-DEPTH RESISTIVITY 
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A) RANGE -0.0008 TO 106n- cm 

B) VOLUME SAMPLED -6x 10-10 cm' 

C) PRECISION ±10% 


G. ROLE OF SUBSTRATE QUALITY ON IC 
PERFORMANCE AND YIELDS 

R. N. Thomas 
Westinghouse R&D Center 
Pittsburgh, PA 


The dramatic progress made in silicon device and IC technology in recent years 
is derived in part from the availability of high quality, large-area silicon substrates. 
The widespread acceptance of the Czochralski growth process by the IC manufac- 
turers arises from the relative ease with which large, doped single crystalline 
material can be produced by crystal pulling from a silicon melt contained in a quartz 
crucible. Over the past decade, engineering efforts have focused almost exclusively 
upon the scaling of crystal growth equipment to larger and larger crystal diameters 
until today; 100-mm diameter crystals have become commonplace. Crystals up to 
150-mm diameter are anticipated in the next few years due to the favorable economics 
of large-area wafer processing (Fig. 11-36). 

Device performance is, however, often limited by the basic substrate properties, 
and there is considerable concern at the present time that the quality of commercially 
available silicon will impose severe limitations as device features are scaled down to 
the micron and sub-micron VLSI range. The adverse effects of crucible dissolution 
on compositional purity combined with micro -segregation effects and crystalline 
imperfections resulting from uncontrolled thermal convection in large volume melts are 
major problems. In silicon LSI, these chemical and structural defects are found to 
interfere with device processing with deleterious consequences to performance and 
yields. 

An even more crucial need for improved materials quality exists in high-speed 
GaAs FETs and the newly emerging monolithic microwave GaAs ICs, which are 
expected to have a broad impact on futxire microwave signal processing and power 
amplification. In the past, the variable and often inferior quality of commercial 
horizontal Bridgeman or gradient freeze substrates has been a major deterrent to the 
emergence of a reproducible, low-cost GaAs manufacturing' technology. Recent 
advances in liquid encapsulated Czochralski equipment offer a capability for producing 
large-area, cassette-compatible GaAs slices of improved quality. 

At the Westinghouse R&D Center, these problems are being addressed through 
crystal growth research of large diameter crystals directed at developing improved 
Si and GaAs by the Czochralski crystal growth process for high performance IC 
processing. 


SILICON 


Infrared Focal Plane Arrays 


Infrared CCD imagers based on ext rinsically- doped silicon form an important 
family of ICs, where performance is strongly influenced by the silicon purity. For 
3 to 5 ym infrared imaging, high photosensitivity depends on indium-doped silicon 
substrates of high indium content and high background purity. Unfortunately, 
dissolution of the quartz crucible in the silicon melt contributes, in addition to large 
amounts of oxygen, boron, and aluminum impurities which are particularly harmful to 
p-extrinsic silicon detectors. 
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The use of high purity synthetic fused silica crucibles in place of the industry- 
standard quartz crucibles in the Czochralski process offers considerable advantage. 
Figure 11-37 illustrates the effectiveness of these synthetic fused silica crucibles. 

An order of magnitude reduction in the boron content is observed. When high purity 
Si:In crystals are compensated by precise neutron transmutation doping techniques, 
extrinsic indium-doped silicon of exceptionally high IR sensitivity is produced (Fig. 
11-38). Undoped silicon crystals produced by high purity Czochralski pulling exhibit 
p-type resistivities in the 3000 to 5000 fi-cra range. 

LSI Silicon 

Oxygen and , to a lesser extent , carbon impurities have a profound influence 
on defect nucleation and precipitate formation in large -area Czochralski silicon. 

Oxygen is present at concentrations in excess of the solid solubility at the tempera- 
tures normally employed in device processing. An oxidation step leads to Si0 2 

precipitation and the generation of oxidation-induced stacking faults (OSF). Carbon 
impurities occur in striated distributions (as do all strongly segregating impurities) 
forming nucleation sites for point defects and interstitial oxygen. These micro-defects 
are sinks to electrically-active impurities introduced inadvertently during processing 
with major consequences to LSI processing. 

Although occurrence of these defects in active junction regions is harmful, 
causing high leakage currents, it is now common practice in complex LSI cycles to 
develop a denuded , low oxygen content surface layer and encourage SiO ^ precipitate 

formation in the bulk for intrinsic gettering (Fig. 11-39). 

Crystal growth studies of large, 100-mm diameter silicon crystals are being 
carried out because of the antid-pated need for substrates of very low microdefect 
content in future very high speed integrated circuit (VHSIC) fabrication. Improved 
purity and uniform incorporation of oxygen are prime goals. The axial and radial 
oxygen distribution in 100-mm diameter silicon crystals pulled in a HAMCO CG200RC 
is shown in Figure 11-40. 

Exploratory studies of magnetic field effects in Czochralski silicon crystal 
growth are also being undertaken. Transverse magnetic fields are found to suppress 
unwanted melt convection (Fig. 11-41) and offer a unique means of controlling the 
oxygen in silicon (Fig. 11-42). 


GALLIUM ARSENIDE 
Monolithic Microwave GaAs Circuits 


The recent emergence of high performance monolithic GaAs circuits for digital 
logic and analog high gain/low noise and power amplification at microwave frequencies 
has caused renewed interest in semi -insulating GaAs crystal growth. 

The use of high purity semi -ins xilating GaAs prepared by LEG growth for 
selective direct ion-implantation processing has been a major factor in the improved 
performance and reproducibility of monolithic GaAs ICs which are currently under- 
going development at different laboratories. Figure 11-43 illustrates round, 50- and 
100-mm diameter GaAs slices cut from <100>-crystals grown in a high pressure liquid 
encapsulated Czochralski (LEC) Melbourn puller. Recent investigations emphasizing 
in-situ synthesis of the GaAs melt in silicon-free pryoltyic boron nitride crucibles to 
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preserve purity, yield reproducible growths of thermally stable, semi-insulting <100>- 
GaAs crystals without resorting to conventional chromium doping. Residual shicon, 

15 - 3 

Group VI , chromium and iron impurities are maintained below the 10 cm range in 
undoped LEG crystals (Fig. 11-44). 


Improvements in the compositional purity of undoped GaAs crystals pulled from 
PEN crucibles are also indicated by the high measured carrier mobility. Mobilities 


of about 6000 cm'^Wsec are observed and correspond to a total ionized impurity con- 


16 -3 

tent in the low 10 cm range (Fig. 11-45). Indirect evidence of the substrate 
purity is also obtained from the transport properties of n-layers formed by Si ion 
implantation of undoped GaAs /PEN substrates. The electron mobility in n-implanted 


1.6 “3 IT “3 

layers with doping levels between 1 x 10^ cm and 5 x 10 cm ^ is found to be 
limited only by impurity scattering and comparison with theoretical predictions again 


indicated about 


1 X 


10 


16 


cm 


residual impurity concentration in the substrate. 


The microscopic inhomogeneity of dopant and residual impurities in GaAs sub- 
strates is cause for concern in large-airea substrate processing of microwave devices 
and monolithic circuits. Impurity striation behavior is observed in large diameter 
LEG GaAs crystals as shown in longitudinal sections of Figure 11-46. In the case of 
semi -insulating GaAs/PEN where the impurity content is low , the observed striations 
are probably the result of local fluctuations in stoichiometry brought about by 
variations in the microscopic growth rate. Measurement of temperature fluctuations 
in a 6-in. diameter, 3 kg E20 2-encapsulated GaAs melt in the high-pressure Melbourn 

system are shown in Figure 11-47. Temperature excursions approaching 10°C are 
observed at the center of the encapsiolated GaAs melt, which are significantly higher 
than those generally observed in unencapsulated Czochralski melts of silicon. 

The improved quality of round, large-area LEG substrates, which are compati- 
ble with modern semiconductor processing equipment , is source of considerable 
optimism that a reproducible, low-breakage GaAs IG processing technology will be 
realized in the near future. 


SUMMARY 

Industrial research and development in silicon and gallium arsenide crystal 
growth are directed towards producing large diameter substrates (because of the 
cost leverage advantages) , which satisfy the specific needs of today’s complex elec- 
tronic devices through the "tailoring" of substrate quality for optimum device proc- 
essing. As the IG technologies are driven more and more towards high operating 
speeds, by scaling of device features to micron and sub-micron geometries, it is 
anticipated that large-area substrates of significantly improved compositional purity, 
dopant distribution and structural perfection on a microscopic as well as macroscopic 
scale will become increasingly important requirements. These needs are currently 
being addressed through the exploratory use of, for example, magnetic fields to 
suppress convection effects in Gzochralski crystal growth. Although the growth of 
large crystals in space appears impractical at present, the current efforts to improve 
substrate quality could obviously benefit from the experiences gained in smaller- 
scale growth experiments conducted in the zero-gravity environment of space. 
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11-36. 
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Figure 11-37. 
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Comparison of responsivity for neutron compensated and 
conventionally compensated Czochralsld Si:In. 
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After Kishino et al, VLSI Laboratories. Japan (1980) 







AXIAL & RADIAL DISTRIBUTION OF OXYGEN 
IN LARGE DIAMETER CZOCHRALSKI SILICON 



Fraction of Melt Solidified 







Effect Of Magnetic Field On 
Temperature Fluctuations In 
Czochralski Si Growth 
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Figure 11-41 



EFFECT OF TRANSVERSE MAGNETIC FIELD ON OXYGEN CONTENT 
IN CZOCHRALSKI SILICON 
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77 


Bulk SIMS analysis of semi-insiilating GaAs prepared by LEG and horizontal 
Bridgeman growth. Note marked reduction of Cr and Si 
in PBN growth runs. (Note scale discontinuity.) 



APPARENT MOBILITY VS TOTAL IMPURITY CONTENT FORN-TYPE 
SEMI-INSULATING GaAs 
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Total Ionized Impurity Concentration, cm 

Figure 11-45. 
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Figure 11-46 






III. POWER NEEDS IN SPACE 
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A. POWER REQUIREMENTS FOR SPACE CRYSTAL GROWTH 


E. Kern 
Consviltant 
Del Mar , CA 


The need for power to grow silicon of a size suitable to test and to eventually 
make into devices must be addressed, so that expected results are tied to power 
which is available. Of particular interest is the projection of size of float- zone 
crystals as a function of time. The accompanying graph (from the book "Floating 
Zone Silicon" by Keller and Muhlbauer) shows that 4 in. (107 mm on graph) will be 
the normally used diameter by 1985 and 5 in. (125 mm) by 1990. Material to be 
tested in device lines in this time frame should be 4 in. or more. By analyzing the 
various power losses, this translated into 25 kW with a 50 percent power efficiency, 
which is much improved over present RF heating efficiencies. Further development 
is therefore needed to provide this improvement - a prerequisite to processing in 
space . 
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1955 1960, 1965 1970 1975 1980 1985 1990 1995 

*-YEAR 

FLOAT -ZONE CRYSTAL DIAMETERS, 1954 TO 1995 

Figure III- 3. 
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B. ELECTRIC POWER NEEDS IN SPACE 


Donald M. Waltz 
TRW 

Redondo Beach, CA 


Mr. Don Waltz of TRW's Defense and Space Group presented the power require- 
ments for some specific float-zone experiments in space. He also gave the present 
available power figures for the Space Shuttle and some projected available power for 
advanced vehicles. The details are shown in the following figures. 
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Figure III- 



FLOAT ZONE PROCESSING ELECTRIC POWER 
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POWER CONVERSION EFFICIENCIES 
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Figure III- 



DATA USED FOR POWER VS. 
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Figure III- 



HOST VEHICLE FLOAT ZONE POWER BUDGETS 
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SAMPLE HEATING ELECTRIC POWER REQUIRED. KILOWATTS 


CADMIUM TELLURIDE POWER REQUIREMENTS 
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CONTROL AND EASE OF EQUIPMENT INTEGRATION INTO THE HOST VEHICLE 



C. THE HISTORICAL TREND IN FLOAT-ZONE CRYSTAL 
DIAMETERS AND POWER REQUIREMENTS FOR 
FLOAT -ZONED SILICON CRYSTALS 


Horst G. Kramer 

Monsanto Industrial Chemicals Co. 
St. Peters, MO 


Mr. Kramer presented an alaysis of the power needed to zone silicon crystals 
by RF heating. The heat loss mechanisms are analyzed in the following figures. 
Curves are then given for power as a function of crystal diameter for commercial 
silicon zoning. 
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Derivation of Equations by A. J. Syllaios and T. N. Lovern II, Monsanto Co, 


Conduction Losses 

The power into the solid over the surface z = o (solid-iiiielt 
interface) is: a 


P , = /“k 
cond / s 


I 8T I 

V 

\ / z = o 


2-rrrdr 


(II-l) 


This is the heat loss from the melt via conduction through the solid. 

Neglecting radial temperature gradients, the one-dimensional steady 
state heat conduction equation is: 




(II-2) 


The heat loss per unit volume due to radiation from the surface is: 


• • • s: — 


o„ 2iradzT 

l 5 S 


2, 

ira dz 


(II-3) 


Thus, 


^ (l. = 

dz \T dz y 


^ „4 

k T a 
m m 


where the expression k = 


k T 

m m has been used for silicon. 


Assuming T = at z = o 


=T=o atz=“ 


Then, 


1c f a“ 
m m 


dT _ ^ ^ 

'^m " y 

T = n 'S ^ 

lir“T~ 

\ ra m > 


h 3 


^B "s 

1c T a 
L m m 


z + T 


-2 -Is 


-2 

Itl 


(II-4a) 


-3/2 

Cll-4b) 


ClI-5) 


The power by conduction is found from eqs. (Il-l)and (Il-4a) 


cond 


= IT 


g e k T ^ ^ a^^^ 

B s m m y 


(II-6) 


Figure III-16. 
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Radiation Losses from the Liquid 


For the vertical surface: 

= 2T. (Og c, i:J) a 


where = zone length, cm 


For the horizontal surface: 




II-7 


II-8 


Contribution of Latent and Specific Heats 
P = 4.18 trdLva^ 

i-l 

P„ = 4.18 TTdC ATva^ 

H p 

where v = zone travel in cm/sec. 

Total Power Requirements 
For disk-shaped zone 

P(disk) = (V) * Pl + 

and for zone formed by needle-eye coil 
P (needle-eye) = P(disk) + 2P^^^ (H) 


II-9 

II-IO 


II-ll 

II-12 


Figure III- 16. Concluded. 



THE VALUES OF THE CONSTANTS USED 


Emissivlty , solid: 

Emissivity, liqixid: 

Melting temperature: 

Thermal conductivity, solid: 

Thermal conductivity, solid, at T^: 
Latent heat of fusion: 

Specific heat , average : 

Stefan- Bolt zman constant: 


e = 0.46 
s 

= 0.3 

T = 1685°K 
m 


= 0.216 w/cm°K 
338 cal/g 

C =0.23 cal/g°K 
P 

- 1 9 9 4 

Og = 5.729 X 10 w/cm^oR^ 


Figure III-17. 
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Figtxre III- 19 





IV. GROUND BASED RESEARCH AND MEA 
EXPERIMENT SUGGESTIONS 




A. SCHEDULE FOR EXPERIMENT PROGRAM DEVELOPMENT 


Iva C. Yates 
NASA -MS FC 
Huntsville, AL 


Mr. Yates from the Materials Processing in Space Project Office presented the 
major phases of the development of facility class experiment hardware. Times typi- 
cally required for both science working group and project activities are shown in 
the following figure. 
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GENERAL SCHEDULE FOR EXPERIMENT PROGRAM DEVELOPMENT 
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Figure IV- 1. 




B. GROUND BASED SILICON ZONING PROGRAM 


E. L. Kern 
Consultant 
Del Mar, CA 


The purpose and goals of the ground-based investigation phase A/B of the 
experimental float- zoning of silicon is given in the following figures. The overall 
program goals, leading to recommending experiments to be done in phase C/D are 
spelled out. Figure IV- 4 lists the program content of the contract with Westech 
Systems, Inc., to start shortly (actual start date was November 17, 1981, and runs 
for 6 months). Figure IV- 7 is a list of thermophysical properties which must be 
accurately known if thermophysical models can be compared to experimental zoning 
of silicon. 


Ill 
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Figure IV- 2. 
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Develop hardware to the point of specifiying flight hardware design 

INCLUDING ALL CONTROLS^ AUTOMATION AND PROGRAMMING AND SHOWING THE 
RELIABILITY OF ■ SUBSYSTEMS UNDER CONDITIONS SIMILAR TO SPACE 
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Figure IV- 
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Figure IV- 7. 




C. PROPOSE MEA EXPERIMENTS FOR ZONING OF SILICON 


E. L. Kern 
Consultant 
Del Mar , CA 


Materials Experimental Apparatus (MEA) experiments are outlined in the follow- 
ing figures. This is a precursory set of experiments which will be proposed to 
NASA, which will result in meaningful investigation at the power levels available in 
early Shuttle flights. 


118 




119 


Figxire IV- 8. 
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Figure IV- 


















D. FLO AT -ZONE EXPERIMENTS 


John Verhoeven 
Iowa State University 
Ames Laboratory 
Ames , OH 


MAIN TASK 

Can surface tension driven convection in float -zone experiments be eliminated 
by means of oxide films (several monolayers)? 


EXPERIMENTAL APPROACH 

1. Select molten Sn as model system because of low vapor pressure and ease of 
handling. 

2. Develop surface analytic techniques to evaluate oxide film thickness on molten tin 
versus temperature, O2 pressure, and time. 

3. Develop a system for float- zone melting of Sn which 

1) Allows surface oxide to be controlled and monitored 

2) Maximizes Marangoni flow and minimizes buoyancy flow 

3) Somehow allows analysis of effect of film on convection. 
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QUALITATIVE DATA ON STRENGTH OF THIN OXIDE FILMS 
W. B. Pietenpol and H. A. Miley (Phys Rev 1927) 


I ^ 1 ->4 cm ► I 


c 


Cu 


COPPER 

ELECTRODE' 


— 0.8 mm DiAM. 
WIRES OF Pb 8i Sn 


Cu 


HEATED IN Oj OR AIR 


I I 



- Strength of oxide sufficient to hold 4 cm of 0.8 mm wire to 600°C 

- Ductility of oxide sufficient to withstand thermal expansion encountered 
by repeated heating and solidification 

- Oxide stable - Held at 250°C above m.p. in air for 1 hour 
Our experience: 

Uranium chunks of 3 mm size with native oxide retain their shape after 

_ 5 

U melts in vacuum of 10 Torr. 
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TEMPERATURE °C 






INTENSITY 



40 45 50 55 60 65 70 75 

ELECTRON ENERGY (eV) 


Energy Loss Spectrum (ELS) of Sn02 



Interband 
Trans i tion 
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Bulk 
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SnO 
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Sn 
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13.5 

9.3 

(3.8) 


Figure IV- 14. 
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mm diam. recess milled 
Into stainless steel plate 



Figure IV- 15. 
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5-Volt peak, AE*4*8tV 
Interbond Oxide peok, AE - 26.3 eV 
Sn02 bulk plosmon, AE - 16 eV 
Sn bulk plosmon, AE ^ 13.5eV 
Interbond Metal peak, AE = 24,25 


■5V(xl/2) 


d liri 



MBP(xl) 



IBM(x1) 


SPUTTER TIME (minutes) 


Figure IV-17. 







CONCLUSIONS 


1. Molten tin surfaces may be readily cleaned by ion bombardment with 2. 4 keV 
argon ions with ion beam current densities of approximately 100 yA7cra2. 

2. It is very important that the ion bombardment beam is positioned with adequate 
clearance to avoid striking the CM A shield. 

3. Molten tin may be held in Ta containers at temperatures up to at least 240°C 
without producing contamination detectable at the tin surface. 

4. The surface of molten tin at 240°C will remain atomically clean for many hours 
at total pressures of 4 x 10“ ^ Torr, and perhaps higher. 
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INCOMING ELECTRON 


I|BM - INTENSITY 



>IBM = jo 


/ 


INTENSITY 
OF IBM SCATTERED 
ELECTRONS WITH NO 
OXIDE (X = o) 




OXIDE THICKNESS X 

AFTER 2 HR OXIDATION = ~ 8 A 


Figure IV- 19. 
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CONCLUSIONS (Concluded) 


6, After holding for around 20 to 40 minutes under vacuum this oxide appears to 
break up into a mobile island structure. 

7. Depth profiling of oxidized molten tin does not appear to lend itself to quantita- 
tive analysis, apparently due to breakup of the film into a mobile island structure. 


IMPLICATIONS FOR OUT STUDY 


1. Will be able to produce thin (10 to 15 ^ continuous oxide films on molten tin at 
O2 pressure low as 5 x 10 " 4 t. 

2. The films will probably only remain continuous if pressure is maintained within 
system to heal cracks. 

3. Only moderate vacuum levels required, ^^^4 x 10"® T. 

4. Will be able to make direct contact to molten zone with Ta shielded T.C. — no 
contamination . 

5. Samples may be cleaned over entire molten surface with a stationary ion beam. 
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DISK 
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Figixre IV-22. 
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EXPERIMENTAL PLAN 


207 

1. Solidify disks with and without oxide films. Measvire solute profiles using Bi 
and look for effects due to oxide. 

If flow velocities are on order of 1 to 10 cmVs and we solidify at rates 
'A^O.OOl cmVs (10 ymVs) we should see an effect. 

2. Look for oscillatory convection on clean zones utilizing miniature thermocouples for 
detection. If find, determine whether oxide will suppress temperature oscillations. 


PROGRESS 

_ g 

1. Check out disk heater in oil pumped (10 T) vacuum system. 

2. Measure temperature profiles and compare to calculations. 
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Initial Work with Disk Heater Design 

- Trouble with Ta/Sn contact 

- Overcome by electropol. Sn disk 



Figure IV-23. 
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RADIAL POSITION (mm) 

Figure IV"24. 


TEMPERATURE 
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Figure IV- 25. 
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THE SHAPE OF INTERFACE 


Figure IV- 26, 
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Figiire IV- 2 9. 
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Figiire IV- 30. 
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ZONE WAMCTEM IN CENTIMETtHE 



MAXIMUM DISK ZONE DIAMETER VERSUS PLATE THICKNESS 


Figiire IV- 32. 
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Figure IV- 33. 
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Figure IV- 34. 
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E. PRESENT GROWTH TECHNOLOGY OF SILICON -GERMANIUM ALLOYS AND 
POSSIBLE ADVANTAGES OF MICRO GRAVITY GROWTH 


Dr. Oscar M. Stafsudd 
UCLA 

Los Angeles, CA 


Dr. Stafsudd, who is Professor of Electrical Engineering at UCLA, reviewed the 
present growth technology of Silicon- germanixjm (Si-Ge) alloys and the possible advan- 
tages of growth in microgravity. Si-Ge alloys are of interest because they have a 
continuous variation of bandgap energy from the germanium bandgap (0.65 eV) to the 
silicon bandgap (1.09 eV) (Fig. IV- 35). The unusual two-slope behavior of Eg versus 
composition is due to the differences in the conduction band structure between Si and 
Ge. Below 17 percent (atomic), the germanium band structure dominates; and above 
it, the bands are "silicon like." 

One can therefore tune the band structure to produce intrinsic photodetectors 
at wave lengths of particular importance, such as 1.06 ym for laser detectors and 
1.3 to 1.6 ym for fiber optics. The competitive materials in the 1.3 to 1.6 ym range 
are 3 to 5 ternary and quaternary compounds, for which the microcircuit technologies 
are not as developed as with silicon materials. Germanium does not possess a native, 
stable, and protective oxide, such as silicon does, thus providing for each processing 
and good electrical protection. It will be interesting to study the oxide system of 
the Si-Ge alloys. 

The extrinsic doping levels of dopant ions in Si-Ge alloys varies with alloy com- 
position (Fig. IV- 36). While this is of interest for infrared detectors, this develop- 
ment has met with limited success. This may partly be due to reduced mobilities and 
other degraded electrical properties, due, in part, to low crystal perfection and poor 
homogenity. 

Silicon and germanium form a complete solid solution, since each have the same 
crystal structures. The phase diagram is shown in Figure IV- 37. The crystal system 
is clearly well behaved, except for the considerable opening of the solidus -liquidis 
curves. This is somewhat misleading, however, since the alloys are extremely diffi- 
cult to grow when 0. 1 < x < 0. 9 in Si Ge^ . 

The problems in crystal growth are similar, in many respects, to several pseudo- 
binary systems of interest, such as HgTe:CdTe. The Si-Ge system is more tractable 
in that it is a true binary system and the components have relatively low vapor 
pressures at the alloy melting point. 

One problem in crystal growth arises from the low (and varying) distribution 
coefficient , as seen for the germanium distribution coefficient as a function of solid 
alloy content (Fig. IV-38). For Si^ 05 ^ 0 ^ half of the germanium in the melt is 

rejected by the solid at the solid-liquid interface. This results in a large variation 
of the germanium content in the liquid as a function of distance in front of the growth 
interface. This presents a serious problem in the growth-interface stability, leading 
to constitutional supercooling. For instance, if one is limited to a temperature gradient 
at the interface of 100°C/in. , the maximum reported stable growth rate is 1/2 mm /hr. 
Uniform instantaneous growth rates in present crystal growers are hard to control in 
this range. 
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The extreme density variation of the alloy with composition leads to further 
complicate the growth. Figure IV- 39 shows the variation of density with alloy com- 
position and Figure IV- 40 shows the variation in lattice parameter with density. If 
one attempts to grow a solid of 50 percent Ge composition, the liquid rejected has a 
density of 4.8 gmVcc, whereas liquid of 50 percent Ge has a density of 3.6 gm7cc. 
This provides a large density drive in this system. 

The effects of this density segregation is seen in Figures IV-41 and IV- 42. 
Figure IV-41 shows the fraction of Ge in a horizontally grown crystal in three vertical 
sections taken along the growth direction. Gravitational separation of the heavier Ge 
leads to a higher concentration at the bottom of each slice. In addition, there is a 
skin of nearly pure Ge over the entire length of the bottom of the crystal. This 
shows one of the few undeniable examples of gravitational separation in melt growth 
at g = 1. 


Figure IV- 42 shows the variation of Ge concentration (vertical scale) as a 
function of growth rate, R. The variation between top and bottom of an ingot is 
greater with faster growth rates. Figure IV- 43 shows how the top to bottom varia- 
tion (vertical scale) as a function of the alloy composition. 

The thermal expansion data for various alloy compositions is shown in Figure 
IV- 44. Si-Ge alloys make efficient thermoelectric generators at elevated temperatures. 
Figure IV-45 shows the thermoelectric figure of merit as a function of composition, 
indicating that the most interesting compositions are also those of high interest to 
detectors (40 percent Ge corresponds to Eg. = 0.88 or 1.2 to 1.3 ym infrared 
radiation) . 

Similar problems in segregation are observed for crystals grown by the 
Czochralski and float -zone growth methods. 

The growth of Si-Ge alloys in microgravity appears very attractive, as pointed 
out by the g ■= 1 growth experiments and the deleterious role of gravity. In particu- 
lar, the float- zone method, in which a liquid zone of controlled starting composition 
(different from the seed and feed rods) , can be used to grow a large amount of useful 
alloy crystal. Large temperature gradients and relatively flat growth interfaces will 
be necessary in order to obtain homogeneous crystals growth. 
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V. DETERMINING THE MARANGONI EFFECT 
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A. REVIEW OF LITERATURE SURFACE TENSION DATA FOR MOLTEN SILICON 


Stephen Hardy 
National Bureau of Standards 
Washington, DC 


Measurements of the surface tension of molten silicon have been attempted by a 
number of workers over the last thirty years (Table V-1). Although the values 

2 

cluster around 730 mJ/m , there is evidence that the surface tension of clean silicon 
surfaces is much higher than this. The work of Reference 2 fonnd that the surface 

2 2 

tension increased from 730 mJ/m to 860 mJ/m in about 15 min in a hydrogen atmos- 

2 

phere suggesting a reduction in oxygen activity. Thus the 730 mJ/m value commonly 

2 

found may lie as much as 130 mJ7m below the true value. 

For Marangoni flow, the important parameter is the variation of surface tension 
with temperature, not the absolute value of the surface tension. There have been 
no measurements of this quantity for silicon. Estimates in the literature are, I 2 

believe, no more than guesses. For instance, a recent book gives d/dT = -0,10 mj/m"^ 
°K and cites References lA and 2 in the table. Presumably this value was calculated 
by taking the difference between the surface tension of Reference lA and that found 
with Zr02 substrates in Reference 2. However, the surface tension value of [lA] 

appears to be a typographical error. The correct value from the original literature 
is listed as [IB]. Note that the use of the correct data would give a positive tem- 
perature coefficient, which is almost certainly wrong for a pure material. It is 
obvious that other pairs of numbers from the table give widely different temperature 
coefficients. In general, it is not possible to calculate temperature coefficients using 
surface tension measurements from different experiments because the systematic errors 
are usually larger than the changes in surface tension because of temperature varia- 
tions. 


The lack of good surface tension data for liquid silicon is probably due to its 
extreme chemical reactivity. Extensive efforts have not found a material which resists 
attack by molten silicon. Thus all of the sessile drip surface tension measurements 
are probably for silicon which is contaminated by the substrate materials. 
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TABLE V-1. 


Temperature 

(°C) 

Atm 

Substrate 

Y mJ/m^ 

Reference 

Method 

1410 

Ar 

- 

730 

lA 

drop weight 

1410 

Ar 

- 

720 

IB 


1450 

He 

ZrO^ 

725 

2 


1450 

He 

Ti02 

730 

2 

sessile drop 

1450 

He 

MgO 

740 

2 


1450 

H 

MgO 

860 

2 


1550 

Vac. 

? 

720 

3 

sessile drop 

1550 

Vac. 

9 

750 

4 

sessile drop 

1550 

Ar 

- 

825 

5 

pendant drop 


lA. Handbook of Physics and Chemistry, Vol. 50, 1969. 

IB. Kek, P. H., Van Horn, W., Phys. Rev. 91, 512, 1953. 

2. Kingery, W. D. and Htimenik, M., J. Phys. Chem., 57, 359, 1953. 

3. Geld and Petrushevski , Izv. A. N. , OTN , 3, 160, 1961. 

4. Dshemilev, Popel, and Zarevski, Fiz. Met. i Met., 18, 83, 1964. 

5. Eljntin, Kostikov and Levin, Izv. Vys Uch Zav. , Tsvetn. Met., 2, 131, 1970. 
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B. CALCULATIONS OF THE SURFACE TENSIONS OF LIQUID METALS 


David G. Stroud 
Department of Physics 
The Ohio State University 
Columbus , OH 


A crucial ingredient in understanding the Marangoni effect and predicting its 
influence in gravity-free environments is a knowledge of the surface tension of liqmd 
metals, particularly as a function of temperature and of impurity concentration. The 
work described in this presentation is aimed at understanding the surface tension of 
liquid metals and alloys from as close to first principles as possible. The two ingre- 
dients which are combined in these calculations are the electron theory of metals, as 
developed by many workers over the last twenty years or so, and the classical theory 
of liquids, as worked out within the framework of statistical mechanics. The results 
of this work are a new theory of surface properties of liquid metals which permits 
one to calculate the surface tensions and surface density profiles from knowledge 
purely of the bulk properties of the coexisting liquid and vapor phases. The method 
works quite well for those pure liquid metals on which it has been tested; work is 
now underway to extend it to mixtures of liqmd metals, interfaces between immiscible 
Hqmd metals, and to the temperature derivative of the surface tension. 

The basis of the theory is described in a recent publication by Mon and Stroud 
[1]. The surface tension is essentially the excess Helmholtz free energy per unit 
area required to form an interface between the liquid and the vapor. Our theory 
divides this free energy into two parts (1) a ’’local" part, arising from the fact that 
in going from the liqmd to the vapor phase the liquid metal must be compressed or 
expanded away from its preferred equilibrixam densities, and (2) a "gradient" term, 
which comes from the fact that a surface forces the liqmd metal to be non-uniform 
when energetically it is favorable for the liquid to be in a uniform state. The actual 
surface profile results from a balance between these two competing terms. We find 
that a high surface tension is correlated with (1) a low compressibility, leading to a 
high resistance to compression or expansion, and (2) a high coefficient of the gradient 
term, which is related in a complicated way to the bonding, between ions in the liquid 
metal. 


We have tested the theory mentioned above by applying it to a group of ten 
simple metals, i.e., metals in which the valence electrons are free-electron-like and 
interact relatively weakly with the ionic cores. The results of this calculation are 
shown in Figures 1 and 2 of this summary, where the surface tensions (experimental 
and theoretical) are plotted against r . r is a linear measure of electron density 

4 3-1 ^ ^ 

and is defined by (^nr ) = n , where n is the average density of valence electrons. 

As may be seen, the theoretical and experimental numbers never differ by more than 
30 percent or so, even though the measxired surface tension differ among themselves 
by more than an order of magnitude. 

Further tests of the method are described in References 1, 2, and 3, where, 
in addition to the pure liquid metals, liquid rare gases are considered. These 
resemble liquid metals in having interatomic forces which are roughly central and 
two-body, and thus are amenable to treatment by the same methods. In all cases, 
the agreement with experiment is quite satisfying, as shown in these papers. 
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Because of this encouraging initial success , we are now in the process of under- 
taking various extensions of this work in the hope of making close contact with 
materials directly relevant to the Materials Processing in Space program. Mon and 
Stroud have extended the theory to treat the interfacial tension between two immis- 
cible liquid metals, and applied the method to treat molten Na-Li alloys. Preliminary 
calculations yield a surface tension of about 10 dyne/cm at a temperature of 500 K, 
where the two metals are somewhat but not perfectly soluble in each other. This 
number is very sensitive to temperature for this particular material. While 10 dyne /cm 
is rather low compared to the surface tensions of the two pure constituents, (2) is 
not unreasonable. 

Further work is underway to calculate the surface tensions of pure liquid metals 
as a function of temperature. Initially, Wood and Stroud are concentrating on liquid 
Na, a metal which is quite well understood in bulk liquid form. It is hoped to under- 
stand precisely what bulk parameters affect not only the surface tension but its 
temperature derivative. If this work is successful, it is intended to consider liquid 
Ga as a next application, with the intention of making contact with the experimental 
work of Hardy and Fine [4], and possible also liquid Si, which is a simple metal in 
its liquid phase and which is potentially of much importance for the gravity- free 
materials-processing program. Also underway is a calculation in which W. H. Shih 
and Stroud are trying to use the alloy form of the theory to understand what kinds 
of impurities will most strongly affect the sxorface tension of various liquid metals, and 
which properties of these impurities are most important in this effect. Initially, this 
work will consider metal-metal mixtures, but it may be possible at some later date to 
extend the work to such important contaminants as oxide overlays or other nonmetals, 
most likely only in a semi-empirical way. 


REFERENCES 


1. K. K. Mon and D. Stroud, Phys. Rev. Lett. 45, 817, 1980. 

2. C. Ebner, W. F. Saam, and D. Stroud, Phys. Rev. A 14, 2264, 1976. 

3. D. Stroud, proposal to National Aeronautics and Space Administration, March 
1981. 

4. S. C. Hardy and J. Fine, NBS report NBSIR 80-202 (July 1980, unpublished). 
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t( dynes /cm) 


Calculated and experimental surface tensions of a number of simple liquid 
metals at their freezing point. The calculated numbers are from Mon and 
Strong, Ref. 1. The results are plotted as a function of r , a linear 

® 3-1 

measure of conduction electron density n . r is defined by (4 tt r 73) 

OS s 



Figure V-1. 
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VI. MODELING AND CHARACTERIZING THE SILICON FLOAT -ZONE 


173 




A. MODELING OF THE FLOAT-ZONE PROCESS 


Robert A. Brown 

Department of Chemical Engineering 
Massachusetts Institute of Technology 
Cambridge , MA 


SUMMARY 


This research is directed towards a fundamental understanding of the inter- 
action of heat, mass, and momentum transport in the floating- zone method for growing 
single crystals from the melt. To accomplish this goal we are developing methods for 
detailed numerical simulation of the transport phenomena in a floating- zone. The 
results of these calciilations can be combined with careful experiments to determine 
the effects of solidification -induced, surface-tension-driven, and buoyancy-driven 
convection in establishing dopant redistribution in the melt and the roles of heat 
transfer in crystal and melt and melt/solid interface shape in determining crystal 
quality. 

To date, state-of-the-art finite -element techniques have been developed for 
efficiently calculating the influence of natural convection in the melt on the shape of 
a melt/erystal interface and dopant segregation in the crystal. These techniques are 
demonstrated here for solidification by the Bridgman technique, another system of 
NASA interest. Extension of these methods to the floating- zone system has begun. 
Numerical techniques have been developed that can calculate the shapes of both the 
meltysohd and melt7gas interfaces simultaneously with the thermal fields in melt and 
solid . 


Models for the fluid flows due to the rotation of the feed and crystal rods have 
been completed and the effects of these flows on dopant segregation are being studied, 
especially in the case of zones longer than can be achieved on Earth. 


MOTIVATION 


1. To provide a theoretical basis for analyzing results of careful experiments where 
limited output is available from the experiment. 

• Macro- and Micro- segregation : Radial and Axial 

• Melt /Solid Interface Shape 
Stirface Temperature (?) 
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MATHEMATICAL MODELLING OF MELT CRYSTAL GROWTH 



PROCESS CONSIDERED 

• Floating-Zone Process: Zero-Gravity Operation 

• Capillary Growth and Zone-Refining of Thin Sheets 

• Directional Solidification 


Figure VI -1. 






REGIMES OF CONVECTION IN MELT GROWTH 
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PICTURE IS VALID FOR: 


• Natural Convection 

• Surface-Tension-Driven Flows 

• Rotation-Driven Flows 


Figiire VI -2 


• STUDIES ARE A COMBINATION OF EXACT ANALYSIS AND DETAILED 


SIMULATION BY FINITE ELEMENT ANALYSIS 


• OUR APPROACH HAS BEEN TO DECOUPLE THE PROBLEM. 
BASED ON 

• Low Prandtl Number 

• High Surface Tension 


-->BOTH VALID FOR SILICON 


STAGES 

1. HEAT TRANSFER ANALYSIS WITHOUT CONVECTION, BUT WITH: 

Melt/Solid Interfaces 
Gas/ Me It Meniscus 

2. CONVECTION AND SEGREGATION 


Figure VI -3. 



THE "COMPLETE" MODEL OF THE FLOATING ZONE 


COMPLICATIONS 


• HEAT TRANSFER * 

• Unknown Shapes of Melt/Solid 
Interfaces 

• Radiation 

• Coupling with Meniscus Shape 
and Crystal Radius 

• CONVECTION 


• Natural Convection* 

• Rotation-Driven Flows* 
« Surface Tension Flows 


Figure VI- 4. 



CONVECTION AND SEGREGATION ANALYSIS FOR DILUTE MELTS 



• EXAMPLE FROM DIRECTIONAL SOLIDIFICATION 


Figure VI -5. 
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PROTOTYPE OF VERTICAL BRIDGMAN GROWTH 
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♦ ADIABATIC ZONE IS DESIGNED TO MAKE INTERFACE PLANAR. 




CONVECTIONLESS GROWTH 
Pr-aOl, A= 1/4 
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TEMPERATURE STREAMLINES CONCENTRATION 
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RADIAL COORDINATE R 




185 


Figure VI-10. 













AXIAL COORDINATE Z 


MELT/SOLID INTERFACE SHAPES 
FOR DESTABILIZED GROWTH SYSTEM 

Pr = 0.01 a =1/4 



RADIAL COORDINATE R 

Figure VI-11. 
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RADIAL SEGREGATION PASSES THROUGH 
A MAXIMUM WITH INCREASING CONVECTION 



Figure VI -12. 



B. FLOAT-ZONE CONTROL SURFACES 


Larry M. Foster 
Science Applications, Inc. 
Huntsville, AL 


The equations of state describing a simplified model of the float- zone process 
are given in Figure VI-13. The process and material parameters in Figure VI-13 are: 


P 


solid Peclet number 


^S 

^S 

M 

T ^ and T 2 

and A2 


T 

h(x) 


liquid Peclet number 

solid conductivity 

liquid conductivity 

solid heat transfer coefficient 

material melting temperature 

solid regions ambient temperatures 

interface parameters, e.g. , the product of the crystal growth 
rate, the liquid density and the latent heat of solidification 

material temperature 

liquid siirface temperature. 


Unfortunately, the abundance of arbitrary boundary conditions in Figure VI- 13 over- 
pose the partial differential equations. 


We propose to replace the boundary condition 



9T 

8r 



(T - Tp 


r - 1 , X < -Q 


by 


T(x,l) = f(x) 


where f(x) is to be computed. The unknown surface control temperature f(x) will be 
constructed so that the resulting set of boundary conditions do not overpose the 
partial differential equations. 
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Figure VI- 13. 
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Figure VI-14. 



PROBLEM: COMPUTE f(x) SUCH THAT THE SOLUTION OF 
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DUAL INTEGRAL TRANSFORM PAIR; 
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REQUIRE: 


( 6 ) 



2 S, 







FOR CONVENIENCE: 


G(x - Q) * E(- x), X < 0 




CONTROL EQ: 
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SOLVE; 



L9' 




NUMERICAL CASE STUDY: (BRIDGMAN-STOCKBARGER) 


Figure VI- 22. 
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VII. SUMMARY 
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SUMMARY OF THE SEPTEMBER 1981 FLOAT-ZONE WORKSHOP 


Edward L. Kern 
Chairman, FZWG 
Del Mar, CA 


Improvements, recommended in the December 1980 Float Zone Workshop, have 
been made in orienting the research of this program toward technologically important 
materials. Silicon, the most important semiconductor crystal, will be the primary 
material iinder investigation. It does create some problems. Because of its high 
melting point and high reactivity with almost all other materials, crucial measurements 
are difficult. Solutions to these problems were discussed and investigators are 
already working on them. 

Industry participation, as recommended in the 1980 Workshop, was the key 
ingredient in developing the rationale for space float -zone experimentation. Device 
related material scientists from several leading U.S. semiconductor companies and 
from U.S. government labs stressed the importance of understanding and making 
improvements in undesired material characteristics which lead to limited performance 
or degradation of devices. The most troublesome, which are largely , unsolved , are: 

1) Resistivity striations 

2) Vacancies, which complex with other impurities. 

The 1980 Workshop had indicated that the direction of the NASA program should 
be to study and control these and radial segregation. Programs which will address 
all of these areas are now beginning. They include eliminating the RF coil slot, 
studying the nature of thermal oscillations, and modeling the heat flows to make the 
interface and isotherms flatter, reducing vacancy formation. Equipment to breadboard 
technology drivers, recommended in the 1980 Workshop, is also being set up. 


SILICON 


The first objective of the float- zone ground and space research program will be 
to understand the growth of silicon single crystals. Silicon will be used as the model 
for all semiconductors. The understanding should be sufficient to lead to growth 
improvements here on Earth (g = 1) , where the bulk of silicon crystals will always be 
grown. 

Why silicon? First, it is the most important semiconductor material to both 
industry and government. Second, there are problems that need to be solved to get 
better important devices such as visible light detectors and imagers , infrared imagers 
and radiation hard space solar cells. Third , one cannot adequately model silicon 
problems with lower temperature crystals. The magnitude of important phenomena 
effecting growth (i.e. , Marangoni driven convection, radiation and buoyancy driven 
convection) are very different at this high growth temperature. 

A secondary payoff, which evolves naturally from the planned program, is the 
growth of uniform crystals hich are needed for standards. The Chairman of the ASTM 
Semiconductor committee indicated the industry /government -wide need for silicon 
measurement standards for: 
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1) 4 point probe resistivity 

2) Spreading resistance 

3) Oxygen determination 

4) Carbon determination. 

Where do we stand on modeling the float -zone system? We have a good, across 
the board, start. Analytical modeling in the axlsymmetric geometry should provide a 
good understanding of the Marangoni convection. Since inaccurate values, such as 
surface tension, can lead to erroneous predictions, a number of thermophysical 
properties may need to be critiqued and remeasured. The entire complex zone is 
being modeled and will be compared with experimentally zoned crystals. For the first 
time, we are asking the analysts to tell us what heat profiles are needed to get a 
desired growth configuration. Usually we empirically develop a heating method, and 
then explain what it gives us. In this case, a lot of nonuniformity! 


SILICON DIAMETERS 


Why do we need large diameter silicon? First, you cannot model a large diameter 
growing crystal (with large molten zone volume) with a small geometry. The modeler 
cannot predict all of the effects accurately. Financial support and the very large 
computer programs have not been provided to model the real shape of a molten zone. 

The second reason for large diameter is to test the improvement that uniform 
silicon has on a. specific device. Modern processing lines and photolithography masks 
are at 2-in. diameter, at the smallest. That will be for proving the material super- 
iority, If devices are to be made in production, 3 in. is standard for float -zone 
crystals and 4 in. will be the "workhorse" from 1983 to 1990. Diameters of Czochralski 
crystals are now 4 and 5 in. 


POWER REQUIREMENTS 


Analyses (3 papers) were in good agreement on the power needed to grow 
silicon at various diameters. These took into account the heat processes in the 
growth system (heat of fusion, specific heat, conduction, and radiation) and the 
power conversion efficiencies of different heating methods. All heating methods were 
reviewed. Two-inch diameter silicon can be fabricated with the power available to a 
pallet experiment in the space shuttle bay (5 to 8 kW), but only if power generation and 
coupling efficiencies are raised from the present 20 to 33 percent efficiencies to 50 
percent (for 8 kW available) or 80 percent (for 5 kW available). Fifty percent 
efficiency will take development effort, but meeting the goal seems probable. Seventy 
to eighty percent would require extensive development. If space produced crystals 
make significant improvements and are wanted for critical devices, then the power for 
3 in. (15 kW at 50 percent efficiency) and 4 in. (25 kW at 52 percent efficiency) 
would need the power of the extended Materials Experiment Carrier (MEC) attached to 
the Space Platform (SP). The 12. 5 kW version is now scheduled for 1988 launch. 

The extended 25 kW version could be in the same time frame as a Silicon MEC unit 
could be developed (1991). 
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SEMICONDUCTOR ALLOYS 


A second semiconductor type which would greatly benefit from microgravity 
processing is semiconductor alloys, such as Si;Ge, InP, As and Hg, CdTe. Si:Ge 
is of interest for military devices, such as designators and imaging systems, in the 
1.0 to 1.5 ym range. Si:Ge and InP, As are of interest in the 1.3 to 1.5 ym fiber 
optics telecommunication area, for the Mgh-band width and low-loss wavelength 
range. These crystals are very hard to grow on Earth because of heavy component 
segregation and constitutional supercooling. Si :Ge, cannot be grown at g = 1 
when 0.1 < X < 0. 9. x i x 

Extremely slow and uniform growth rates are needed in microgravity. These 
are also conditions to grow Si:X, where X is a heavy dopant (Tl) needed in high 
concentrations for higher operating temperature infrared imagers. These are of 
interest to the Air Force. 


PRECURSORY EXPERIMENTS 


The rationale for experimentally growing silicon in space, at larger diameters ^ 
is being developed. Is there a reason to do early experiments in space? Yes. In 
a complex growth system, one cannot always predict what wiU take place. For 
instance, the cooling or insulating effect of an inert zoning gas (argon) cannot be 
predicted well, and slow growth rates in micro gravity for Si:Ge cannot be modeled. 

Possible Materials Experiment Apparatus (MEA) experiments were proposed and 
the scientific background surveyed. The workshop reviewed and discussed these 
areas. The doped silicon and Si:Ge alloys will be proposed shortly. This would 
provide an early trial for testing conditions for zoning silicon. 

Discussion led to a suggested experiment for a self-contained cannister con- 
figuration, probably on a sounding rocket or a Get -Away- Special. It involves a 
simple experiment of melting the center of a silicon slice and observing (1) the shape 
of the solid -gas -liquid trijunction and (2) the magnitude of Marangoni flows. 


CONCLUSION 


It is now probable that a well-planned ground based investigation can proAude 
the needed modeling to predict the effect process changes will have on crystal growth, 
including the change of gravity. Modeling and corresponding experimentation can be 
complemented by precursory experiments in space. This program should provide the 
rationale for proceeding with a major space experiment. Definition of a major space 
experiment should be possible in 1984. 
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